Determination of magnetic order of the rare-earth ions in multiferroic TbMn2 05 
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We have employed resonant x-ray magnetic scattering to specifically probe the magnetic order of 
the rare-earth ions in multiferroic TbMn2 05. Two energy resonances were observed, one originated 
from the El-El dipolar transition and the other from the E2-E2 quadrupolar transition. These 
resonances directly probe the valence 5d band and the partially occupied 4f band, respectively. 
First, full polarization analysis, which is a measurement of the scattered polarization as a function 
of incident polarization, confirmed a spin polarization of the terbium valence states (probed by the 
El-El transition) by the Mn^^ spin density in the commensurate phase. Second, full polarization 
analysis data were collected in the low-temperature incommensurate and commensurate phases when 
tuned to the E2-E2 resonance. By employing a least-squares fitting procedure, the spin orientations 
of the terbium ion sublattice were refined. 

PACS numbers: 78.70.Ck, 75.30.Gw, 75.50.Ee, 75.47.Lx 



I. INTRODUCTION 



A small group of materials known as magnetoelectric 
multiferroics has been of major scientific interest for some 
yearsii^i^i^. These materials display phenomena in which 
both magnetic and electric orders are coupled within a 
single phase. In particular one of the most dramatic ef- 
fects has been observed in TbMn205, a system in which a 
complete reversal of electric polarization is induced in ap- 
plied fields of approximately 2 T.^ The mechanism driv- 
ing the multiferroicity in such systems is yet to be fully 
understood; although in studies of the i?Mn205 series {R 
= rare earth, Y or Bi)^ it is clear that the rare-earth ions 
play an important role in the magnetoelectric coupling. 
It has recently been shown by Koo et al& that numerous 
magnetic orders exist in the low-temperature phases of 
TbMn2 05, centered on Mn^+, Mn^+, and Tb^+ sites, all 
with the same wave vector, giving rise to a complex spin 
configuration. The terbium magnetic order couples with 
the lattice, the effect of which increases at low tempera- 
ture due to an increase in the magnitude of the average 
terbium magnetization. It is however not possible to di- 
rectly probe the magnetic moment of the rare-earth ion 
using neutron-diffraction or bulk transport or magneti- 
zation measurements. 

In the past two decades, magnetic x-ray scattering has 
proved to be a useful tool in measuring long-range mag- 
netic order in single crystals. Upon tuning incident x- 
rays to an absorption edge, the scattered intensity from 
the very weak magnetic signal can be dramatically en- 
hanced. This, coupled with the brightness of third gen- 
eration x-ray sources, makes it possible to observe other- 
wise undetectable weak reflections from long-range mag- 
netic, charge and orbital oideiii^^^^ In the case of 



magnetic scattering at the terbium L edges, the signal is 
~ 10^ times stronger than found off resonance. In this 
paper we employ full polarization analysis to investigate 
the polarization dependence of the scattering amplitude, 
with the aim of refining the direction of magnetic mo- 
ments in specific terbium electronic states in TbMn205. 

Full polarization analysis is a relatively new technique, 
which has been shown to be capable of unravelling mul- 
tipole resonances in K2Cr04^^ and modeling competing 
magnetic domain contributions within NpRhGas.^^ The 
Appendix briefly outlines the theory behind the scatter- 
ing amplitude's dependence on magnetic moment vector 
and polarization, providing the basis for the simulation 
and model refinement of the polarization dependences 
measured by us. 

TbMn205 crystallizes into the space group Pbam with 
lattice parameters a = 7.3251 A, 6 = 8.5168 A, and 
c = 5.6750 A, as measured by neutron scattering— li^ii^. 
Mn^+ and Mn^+ ions sit in octahedral and square based 
pyramid oxygen coordinations, respectively^^. 

At low temperature the manganese sublattice is found 
to exist in incommensurate (ICM) and commensurate 
(CM) magnetic phases, in which for both cases the man- 
ganese ion moments align in the ab plane forming two 
spin density waves with wave vector ((5, 0, r), reliev- 
ing an otherwise geometrically frustrated systenJ^ii^. 
Along the c-axis these moments alternate with ferromag- 
netic and antiferromagnetic layers^^. At Tn = 43 K, 
the system enters an incommensurate phase (ICM2) in 
which manganese ions order antiferromagnetically with 
S 0.5 and r = 0.3^^. Slightly lower in temperature, 
Tfe = 38 K marks the onset of ferroelectric order with 
P II b^. The symmetry of the system in this ferroelectric 
phase is reduced, probably to space group Pb2im^'^^ . In 
the temperature range 33 K > T > 24 K the system locks 
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FIG. 1: (Color online) Illustration of the low temperature 
magnetic structure of TbMn205 in the a6-plane, adapted 
from Fig. 3, Chapon et al:^. Antiferromagnetically aligned 
moments lie on the Mn*^^ (pink) and Mn^^ (green) ions that 
are located in square based pyramid and octahedral coordi- 
nations respectively; the terbium ions are shown in purple. 
Oxygen ions are omitted for clarity. 



into a CM phase, (^^ = ^, t = |), which is diagrammati- 
cally shown in Fig. (TJ as refined from neutron diffraction 
experiments. ^^^^^ At T < 24 K TbMn2 05 enters another 
incommensurate phase (ICMl), {S = 0.48, r = 0.32)^^ 
an unlikely transition as one might expect a commensu- 
rate ground state of any simple ordered system. 

Below T/v, the manganese magnetic structure spin po- 
larizes the electronic states of terbium ions. A phe- 
nomena shown to occur and discussed not only in 
TbMn2 05^^, but also in other manganese oxide multi- 
ferroics such as TbMnO^ and HoMn203i^. Blake et 
al}^ suggested that this is specifically due to an inter- 
action with the Mn^+ spin density. The existence of 
a separate magnetic order on the terbium ions was hy- 
pothesized with the terbium sublattice ordering below 
10 K.^^ In this paper we investigate the interaction with 
the manganese spins and any independent ordering of the 
terbium ion sublattice by resonant x-ray magnetic scat- 
tering (RXMS) full polarization analysisi^iii^. It can be 
assumed that the terbium 5d valence band has a large 
overlap not only with the lower energy terbium 4f band, 
but also with the 3d manganese band. The terbium 5d 
states will therefore be spin polarized via interaction with 
neighboring Mn^+ , Mn^+ , and Tb electronic spin density. 
At the terbium Lni edge, the 5d band is probed by the 
2p o 5d dipole transition. The 4f band of terbium, in 
which the unpaired electrons exist, is probed at this edge 
by the 2p <^ 4f quadrupole transition. 



II. EXPERIMENT 

A high quality single crystal of TbMn205 with dimen- 
sions approximately 2 mm was grown at the Department 
of Chemistry of the National Taiwan University by fiux 
growth. The single crystal was prepared such that the 
(4, 4, 0) Bragg refiection was close to the surface normal. 
The sample was mounted on the six-circle diffractometer 
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FIG. 2: (Color online) The transmitted intensity in the cr'- 
channel as a function of the phase plate theta angle. The 
diamond crystal behaves as a half-wave plate when aligned 
such that theta corresponds to a minimum in the graph as 
shown by the blue dashed line^°. This is symmetrical about 
the Bragg angle. 



in EH2 on beamline ID20 at the European Synchrotron 
Radiation Facility with the c-axis perpendicular to the 
horizontal scattering plane. A helium fiux cryostat was 
used to achieve a stable sample temperature through the 
different phases. 

The technique of using a phase plate, positioned in 
the incident beam such that it rotates the horizontally 
polarized light to any arbitrary linear polarization, was 
employed. A diamond crystal was used as the phase re- 
tarder by scattering near the (111) refiection. A conve- 
nient forward-scattering geometry is possible due to the 
low absorption coefiicient of diamond. The phase plate 
was calibrated by measuring the scattered intensity in the 
a channel when rotating the phase plate theta through 
the Bragg angle, see Fig. [2l A minimum in the cali- 
bration curve (blue dashed line in Fig. [2]) was selected as 
the angle at which the diamond crystal behaves as a half- 
wave plate. Note that a x phase plate rotation about the 
incident beam (see Fig. [3]) equates to a e = 2% rotation in 
polarization. The scattered beam polarization was deter- 
mined using an Au(222) analyzer crystal. This technique 
of full polarization analysis is described in greater detail 
in the Appendix and is diagrammatically illustrated in 
Fig. 13 

The polarization of the scattered beam was measured 
as a function of incident polarization. By simulating the 
response and comparing it with the data, it was possible 
to refine the magnetic moment direction. 



III. RESULTS AND DISCUSSION 

The first harmonic magnetic satellite refiection (4+(5, 
4, 0-r) was observed in the ICMl phase (Fig. HJ inset). 
All the scattered intensity was found to be in the rotated 
TT — a' channel, as expected of the magnetic scattering 
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FIG. 3: (Color online) A schematic of the experimental setup. 
The polarized x-ray beam (blue) enters from the left. The 
polarization of the incident beam was rotated by a diamond 
phase-plate. The scattered beam polarization was determined 
through a secondary scattering process from an analyzer crys- 
tal prior to detection. Note the a and tt orientations where 
primes (eg. a\ tt') refer to the scattered beam. 



cross section. The peak, measured along the [100] direc- 
tion, shows a good fit with a Lorentzian squared function 
centered at /i = 4.483. The inverse correlation length, de- 
fined as: 



a 



(1) 



where a is the lattice parameter and Ah is the half width 
at half maximum, gave = 2.212 x 10~^ ± 6.2 x 
10-^A-i. The inverse correlation length of this reflec- 
tion remained constant within experimental error within 
the temperature range 2<r<35 K, except for a slight 
increase close to the ICMl/CM transition. 

Scans of the scattered intensity as a function of in- 
cident x-ray energy through the Tb Lni edge at con- 
stant wave vector were performed at the (4+^, 4, 0-r) 
reflection in both the incommensurate and commensu- 
rate phases. Two distinct energy resonances were ob- 
served, measured in the i\ — a' channel in both the low- 
temperature incommensurate and commensurate phases, 
and are shown in Fig. 2^ and Fig. [Dd, respectively. These 
excitations, which are evident in both phases, are cen- 
tered at 7.510 and 7.518 keV. In comparison with the 
fluorescence spectrum as shown in Fig. IIJ), the higher- 
energy peak occurs just above the Tb Lni absorption 
edge, as one would expect for an El-El dipole transi- 
tion. The other resonance is likely to be of quadrupolar 
E2-E2 origin. It exists 8 eV lower, a characteristic shift 
between dipole and quadrupole energy resonances in the 
rare-earth series^^^^^^^^. This resonance has been studied 
in magnetic field by Bland et al?^ 

The wave vector positions in reciprocal space and the 
transition temperatures of the ICMl and CM reflections 
were measured and compared with neutron diffraction 
results to provide a good indication as to whether or not 
the reflections are of magnetic origin. Figure [5] shows the 
position in h [Fig. [Sfa)] and position in / [Fig. [Sfb)] as 
functions of temperature for the (4+^, 4, 0-r) peak. The 
transition from an incommensurate to a commensurate 
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FIG. 4: a) A scan of the scattered intensity of the (4.48, 4, - 
0.32) reflection as a function of the incident x-ray energy mea- 
sured at constant wavevector at 2 K within the ICMl phase, 
measured through the Tb Lm edge in the tv — a' polarization 
channel. The inset shows the peak profile as measured along 
h. b) A similar scan of the (4.5, 4, -0.25) reflection measured 
at 25 K within the CM phase, c) The measured fluorescence 
curve of TbMn205 showing the Tb Lm absorption edge at 
7.516 keV. 



wave vector on warming is clear and agrees well with pub- 
lished 5 and r values In addition we observed that 
the two phases coexist in the temperature range of 20 - 
22 K. This coexistence has also been observed by both 
soft x-ray and neutron diffraction studies^^^^^, showing 
the phase transition to be strongly first order. The in- 
verse correlation length as measured by the width of the 
reflection increases in the coexistence region, evidence 
of greater disorder of the terbium magnetic structure. 
Figure [Sb shows the temperature dependence of the inte- 
grated intensity of the (4+(5, 4, 0-r) reflection measured 
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Figure [6] shows the Poincare- Stokes parameters PI and 
P2, which are defined as 
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FIG. 5: (Color online) The wavevector (position in reciprocal 
space) measured along h (a) and / (b) of the (4+5, 4, O-r) 
magnetic satellite reflection as a function of temperature in 
the range 2 K < T < 34 K. The CM - ICMl phase transi- 
tion at approximately 21 K is clearly seen, together with the 
phase coexistence occuring between 20 - 22 K. Within the 
ICMl phase the incommensurate wavevector varies in both 
h and / as a function of temperature, c) The temperature 
dependence of the integrated intensity of the (4+5, 4, O-r) 
reflection measured in the h direction. The black (ICMl) and 
blue (CM) data were collected on different experiments. Di- 
rect comparison of the scattered intensity in the two phases is 
therefore not possible, however one can observe the trend of 
the CM reflection decreasing in intensity on decreasing tem- 
perature, and the growth of the ICMl reflection on cooling 
below 22 K. 



in the h direction. Upon cooling, the intensity of the CM 
reflection decreases. At 22 K, the ICMl reflection simul- 
taneously begins to increase in intensity and continues to 
do so below 20 K, at which temperature the CM reflec- 
tion becomes extinct. The CM and ICMl data, which are 
shown in the flgure in blue and black, respectively, were 
measured on different experiments. It is therefore not 
possible to draw direct comparisons of the magnitude of 
the integrated intensity in the two phases. We note that 
the widely acknowledged transition temperature of 24 K 
falls outside this range. This discrepancy is likely to be 
a result of x-ray beam heating or a discrepancy in the 
experimental calibration due to the limited proximity of 
the temperature sensor to the sample. 

Full polarization analysis of the (4+5, 4, O-r) reflection 
was performed at an energy of 7.518 keV corresponding 
to the El-El transition, probing the terbium 5d band. 
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measured as functions of incident polarization in the CM 
phase. It was hypothesized that the terbium ions are 
polarized by the manganese magnetic structure, in par- 
ticular by the Mn^+ spin density^^. As explained in the 
introduction, by probing the terbium 5d valence states 
at this transition, any reflnement or simulation of mag- 
netic structure would reveal the origin of the polarization. 
Superimposed in Fig. [6] are three simulations, the blue 
solid line shows the line shape expected if the terbium ion 
magnetism is solely a result of interaction with the Mn'^^ 
spin density and the dashed red line is the line shape ex- 
pected if the scattering occurs due to interactions with 
the Mn^+ spin density. The green dash-dotted line is 
the lineshape of a simulation of scattering due to inter- 
action with terbium 4f magnetic moments, the directions 
of which were reflned by the flt to the data taken at the 
E2-E2 transition presented later in this paper. This Tb 
4f magnetic structure also exists as a consequence of spin 
polarization by the manganese magnetic structure and 
will itself interact with the 5d band. The reduced val- 
ues of the Mn'^+, Mn^+, and Tb simulations are 2.8, 17.9, 
and 5.8, respectively. This result conclusively verifles the 
hypothesis of Blake et al}^ that the terbium sublattice 
is polarized by the close proximity of the Mn^+ spin den- 
sity, as opposed to Mn^+ spin density. A reduced 
value of 2.8 of the Mn^+ simulation, when compared to 
the value of 5.8 of the Tb 4f band simulation, also shows 
that the Tb 5d band is predominantly polarized by the 
Mn^+ 3d band; however, one clearly cannot exclude the 
interaction with the terbium 4f spin conflguration. In- 
deed, in a RXMS study of HoMn205^^ performed at the 
dipole transition at the Ho Lni edge, an azimuthal depen- 
dence on scattered intensity measured in the commensu- 
rate phase, showed excellent agreement with a theoretical 
azimuthal dependence on scattering from holmium with 
the 5d band polarized by the 4f band magnetic structure. 

To conflrm this result, a conventional azimuthal de- 
pendence of the scattered intensity in the cr — it' channel 
was measured at the same energy at the (4+(5, 4, 0+r) 
reflection in a vertical scattering geometry. The data, 
shown in Fig. [7J compared to a simulation of the ex- 
pected azimuthal dependence, assuming scattering from 
the terbium sublattice when polarized by Mn^+ spin den- 
sity (solid line in Fig. [7j), also shows excellent agreement. 

The polarization analysis was performed again at the 
same energy (7.518 keV) in the incommensurate phase. 
The data is shown in Fig. [HI Due to a lack of a magnetic 
structure reflnement in the literature of the ICMl phase 
the large number of free parameters results in an unre- 
liable reflnement using this technique. We can therefore 
only make comparison to simulations of known magnetic 
structures. The CM Mn^+ magnetic structure simula- 
tion, which is shown to be of best agreement to the CM 
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FIG. 6: (Color online) A plot of the measured Poincare- Stokes 
parameters PI and P2 as a function of the incident x-ray 
polarization of the commensurate (4+^, 4, 0-r) reflection at 
the El-El energy resonance at 25 K. Simulations of the Mn^^ 
and Mn^^ magnetic structures [as refined by Blake et al^^] 
are shown as blue solid and red dashed lines respectively. A 
simulation of the terbium magnetic structure scattering at the 
El-El transition in this phase, as refined later in this paper 
at the E2-E2 transition, is shown by the green dash-dotted 
line. 



polarization analysis presented in Fig. [6l is shown here 
as a blue dashed line. In comparison to the data, the 
sensitivity of this technique to a magnetic structure rear- 
rangement upon the CM/ICMl transition, resulting in a 
change of spin polarization of the 5d band, is clear. This 
could be due to a realignment of the manganese magnetic 
structure or a stronger interaction with the terbium sub- 
lattice or both. The green dash-dotted line plotted in 
Fig. [8] is a simulation of the scattering expected due to 
interaction of the Tb 5d band with the terbium mag- 
netic structure, as refined later in this paper at the E2- 
E2 transition. From the poor agreement of the terbium 
simulation, we would conclude that in this phase the 5d 
band remains predominantly polarized by the manganese 
magnetic structure. The difference in full polarization 
analysis data between the two phases is therefore likely 
be a result of a realignment of Mn^+ moments. 

The polarization analysis measurement was conducted 
in the commensurate phase at the E2-E2 resonance 
(7.510 keV). At this transition we directly probe unpaired 
electrons in the terbium 4f band and, hence, the ter- 
bium sublattice magnetic structure. Figure [9] shows the 
data and a least-squares fit (blue solid line), similar to 
the analysis used in a study of UPds^^. There are four 
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FIG. 7: The measured intensity of the commensurate (4+^, 4, 
0+r) reflection at the El-El energy resonance at 25 K in the 
a — 7i' channel as a function of azimuth angle of the sample. 
The solid line shows a simulation of the azimuthal dependence 
expected when the terbium magnetic ions are influenced solely 
by the Mn'^^ magnetic structure [as refined by Blake et al.-^]. 



crystallographically distinct terbium ions in the magnetic 
unit cell, labeled here according to the convention used 
by Blake et a/.,^^ as shown in Fig. [TTl In a recent 
neutron single crystal diffraction study of HoMn205^^, 
which shows very similar macroscopic properties and mi- 
croscopic structure to TbMn205, holmium ion positions 
1 and 2 were found to have the same moment direction 
as were holmium ions 3 and 4. This is in contradiction 
to the earlier neutron powder diffraction study by Blake 
et al}^ in which both Tb and Ho compounds have the 
rare-earth positions 1-4 and 2-3 paired in the refinement. 
In fitting the data it was evident that it was necessary to 
pair the terbium ions as in the recent HoMn2 05 single 
crystal study. Fixing the a component of moments on 
sites 1 and 2 and the b component of moments on sites 
3 and 4 gave a fit with four degrees of freedom. Mag- 
netic moment directions on the terbium 1 and 2 sites 
were refined to be 10.5 ± 2.6° in the a6-plane relative to 
the a-axis and 0.2^0. lo out of plane, and on the terbium 
3 and 4 sites 292. 5^2. Oo in the a6-plane relative to the 
a-axis and 0.2^0. lo out of plane. 

These moment directions are in contradiction with 
neutron powder diffraction measurements on the 
i?Mn205 series as reported by Blake et a/.,^^ however, 
so are recent refinements of HoMn205 by neutron and 
synchrotron single-crystal data^^^^^^^^ when compared to 
the same 2005 study. In these HoMn2 05 experiments, 
holmium directions were reported to be pointing along 
the a and b axes with a small component in the direction 
of the c-axis. These directions are comparable to those 
refined in this study on TbMn2 05. Here, terbium mo- 
ments are found to be coupled in the same fashion and lie 
approximately along the crystallographic axes in the ab- 
plane (the moments are in fact refined to be canted away 
from the axes by approximately 10°) with a slight c-axis 
component. However, the configuration of the spin ori- 
entations of the rare earth and the Mn^+ magnetic struc- 



6 




20 40 60 80 100 120 140 160 180 20 40 60 80 100 120 140 160 180 

Incident Polarization [Degrees relative to jt ] Incident Polarization [Degrees relative to jt ] 



FIG. 8: (Color online) A plot of the measured Poincare- Stokes 
parameters PI and P2 as a function of incident x-ray polariza- 
tion of the incommensurate (4+^, 4, O-r) reflection at the El- 
El energy resonance at 2 K. The commensurate Mn^^ mag- 
netic structure simulation shown in blue in Fig. [6] is shown 
here as a blue dashed line. The green dash-dotted line shows 
a simulation of the terbium magnetic structure scattering at 
the El-El transition in the ICMl phase, as refined later in 
this paper at the E2-E2 transition. 



ture differ to that of the holmium compound. As a conse- 
quence, it is not possible to explain the terbium ordering 
in terms of a simple antiferromagnetic superexchange in- 
teraction between terbium and nearest-neighbor Mn^+ 
ions where the manganese moment a and h axis compo- 
nents cancel at sites 1-2 and 3-4, respectively, the argu- 
ment which was put forward by C. Vecchini et al?^ in 
the HoMn205 study. 

The ICMl phase of the Rybi20^ series is less under- 
stood, with no complete neutron or synchrotron study 
of the magnetic structure published to date. Full po- 
larization analysis performed below 10 K at the E2-E2 
transition will refine the moment directions of the ter- 
bium ion magnetic order. ^ These can then be compared 
to those refined in the commensurate phase. 

Figure [10] shows the Poincare- Stokes parameters PI 
and P2 measured as a function of incident x-ray polar- 
ization at 2 K in the ICMl phase. Superimposed is a 
least-squares fit (red solid line), assuming the same pair- 
ing of terbium as in the CM phase. Again the a and h 
components of moments on sites 1-2 and 3-4, respectively, 
were fixed, giving a fit with four degrees of freedom. The 
refined moment directions on the terbium 1 and 2 sites 
were 322. 0^4. 6o in the a6-plane relative to the a-axis and 
0.5^0. lo out of plane, and on the terbium 3 and 4 sites 



FIG. 9: (Color online) A plot of the measured Poincare- Stokes 
parameters PI and P2 as a function of incident x-ray polar- 
ization of the commensurate (4+5, 4, O-r) reflection at the 
E2-E2 energy resonance at 25 K. A least squares fit with a 
reduced value of 13.9 is shown in blue (solid line). See 
text for the refined moment directions. The incommensurate 
fit performed in Fig. [10] is shown by the red dashed line. 

it was 308.9^ 10. 4o in the a6-plane relative to the a-axis 
and 0.6^0. lo out of plane. This refinement of the 2 K ter- 
bium sublattice magnetic structure results from Tb-Tb 
ordering below 10 K as predicted^, as well as interac- 
tion with the Mn^+ sublattice, the magnetic structure of 
which is unknown in the ICMl phase. However below 10 
K, the terbium order dominates the bulk magnetization 
of the sample^. 

Also plotted as dashed lines in Fig. [9] (CM) and Fig. 
[lO] (ICMl) are the ICMl and CM fits, respectively. On 
visual inspection there appears to be little difference be- 
tween the two fits, particularly when comparing them 
with the data in Fig. [TOl It is surprising then to find that 
the refinement suggests significantly different spin config- 
urations. The fitting process is clearly sensitive to small 
changes in the the polarization dependence and there- 
fore demands high experimental accuracy. This poses an 
important limitation of the technique, one that must be 
considered when interpreting the results. However, due 
to the small error on the data, we have refined the mo- 
ment directions with a high degree of precision. 

The Mn^+ and Mn^+ magnetic structures have been 
shown not to change on substitution of different E?^ 
rare-earth ions in the commensurate phase.^ This 
strongly suggests that the manganese ion magnetic struc- 
tures are responsible for the magnetoelectric coupling 
in the i?Mn205 series. Experimental confirmation that 
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FIG. 10: (Color online) A plot of the measured Poincare- 
Stokes parameters PI and P2 as a function of incident x-ray 
polarization of the incommensurate (4+^, 4, O-r) reflection at 
the E2-E2 energy resonance at 2 K. A least squares fit with 
a reduced value of 7.5 is shown in red (solid line). See 
text for the refined moment directions. The commensurate 
fit performed in Fig. [9] is shown by the blue dashed line. 



the terbium ions are spin polarized specifically by the 
Mn'^^ spin density is therefore of great significance in 
further understanding of the underlying coupling mecha- 
nism, and, by extension, the magneto-electric properties 
of these compounds. 

Probing the terbium ion 4f band in the CM and ICMl 
phases at 25 and 2 K, respectively, has provided a direct 
measurement of the terbium sublattice spin configuration 
in both phases. Importantly, at 2 K we have observed the 
terbium sublattice ordering itself. By assuming a pair- 
ing of terbium ions as found by neutron single-crystal 
data on HoMn2 05,^^ precise moment directions of the 
probable terbium magnetic structures were found. We 
observed a change in spin configuration due to a weaken- 
ing of the terbium ion average magnetization as expected 
above 10 K,^ allowing for a stronger interaction with the 
Mn^+ spin density, as shown by the dipolar resonance 
that probes the 5d valence band. Figure [TTl diagrammat- 
ically illustrates the terbium moment directions refined 
in both the CM [Fig. [nja)] and the ICMl [Fig. [njb)] 
phases. The moment directions on manganese ions have 
been omitted in Fig. [TTb as they are unknown in the in- 
commensurate phase. By visual comparison of the polar- 
ization dependence we observe little change in CM/ICMl 
transition; however, when fitting to the data we refine a 
significant change in spin configuration with precise er- 
ror bounds. As previously discussed, this illustrates a 




FIG. 11: (Color online) Diagrammatic illustration of the re- 
fined terbium ion magnetic moment directions (black arrows) 
in a) the CM phase and b) the ICMl phase. There exists 
a slight component out of the ah plane not shown here. The 
manganese ion magnetic moment directions in the commensu- 
rate phase are shown by grey arrows as in figure [T] These di- 
rections are currently unknown in the incommensurate phase. 

key limitation of the technique. 



IV. CONCLUSIONS 

Two energy resonances have been experimentally ob- 
served at the (4+(5, 4, O-r) reflection, one at the Tb 
Liii edge while the other at 8 eV lower, in both the 
low-temperature (ICMl) and (CM) phases of TbMn205. 
The higher energy resonance originates from the El-El 
dipolar transition, which probes the terbium 5d band, 
and the lower energy resonance originates from the E2- 
E2 quadrupolar transition, which probes the terbium 4f 
band. The temperature dependence and wave vector 
measurements confirmed the reflection to be of magnetic 
origin. Full polarization analysis conducted at the El-El 
resonance confirmed that the 5d band of terbium is polar- 
ized predominantly by the Mn^+ spin density in the com- 
mensurate phase, with an additional interaction with the 
spin polarized Tb 4f magnetic structure. Tuning to the 
E2-E2 excitation at 2 K (ICMl) and fitting to full polar- 
ization analysis data enabled the refinement of directions 
of the spin ordering of the unpaired terbium 4f electrons. 
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The measurement and refinement were repeated in the 
CM phase. 
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where k is the x-ray wave vector, e is the x-ray linear 
polarization orientation, and z is the unit vector in the 
direction of the ion's magnetic moment, ^ei-ei 

are coefficients dependent upon the transition 



,(1) 

E2-E2 



in question, as defined by Hannon et al.^'^ and Hill and 
McMorrow^^, that determine the strength of the reso- 
nance. Primes refer to the scattered beam. The polar- 
ization, magnetic moment and scattering amplitudes are 
defined using the common coordinate system as defined 
by Blume and Gibbs^^(see Eqs. IA.7p . 

The total scattering amplitude incorporates a phase 
factor and is written as 



RXMSJci.r, 

E1,E2 ^ 



(A.3) 



APPENDIX: FULL POLARIZATION ANALYSIS 

The scattering amplitude, f^^^^ ^ of resonant x-ray 
magnetic scattering is anisotropic; dependent upon the 
direction of magnetic moment. It contains components 
from both the dipole, El-El (see Eq. lA.ip . and the typi- 
cally less intense quadrupole, E2-E2 (see Eq. IA.2p . tran- 
sitions. The E2-E2 transition usually occurs a few eV 
below the absorption edge, due to a stronger interaction 
between the core hole and excited electron. The struc- 
ture factors are extensively derived from spherical har- 
monics by Blume and Gibbs,^^ Hannon et a/.^^, and Hill 
and McMorrow^^. The theory presented here is minimal 
compared to that described in the aforementioned origi- 
nal literature. Only the final results useful in data analy- 
sis are given, that is, the scattering amplitude expressed 
in terms of polarization and magnetic moment direction. 
The formulation of the following equation o^^i^^ (Eqs. lA.ll 
and IA.2P is based on the assumption that the system is 
isotropic where only the electron spin breaks the sym- 
metry. In the case of magnetically frustrated multifer- 
roics, of which TbMn205 is an example, this assumption 
is an over simplification. However the scattering ampli- 
tude in the form presented here has been successfully 
employed, in particular, to the closely related compound 
TbMnOa^^. There are three terms that contribute to the 
El-El resonant scattering amplitude and five terms that 
contribute to the E2-E2. In this experiment, we mea- 
sured first harmonic magnetic satellite reflections and are 
therefore only concerned with the first order terms i^^i^^ 
At the El-El transition. 



rRXMS _ 
JEl-El — ^-^ El-El 

and at the E2-E2 transition, 



(1) 



(e' X e) 



RXMS 
E2-E2 



iF^y_^2[(k'-k)(e'xe)-z 
(k' xk)(e'-e)-z] 



(A.l) 



(A.2) 



where q is the scattering vector and r is the crystallo- 
graphic coordinate of the jth terbium ion. 

The incident x-ray polarization was rotated by a di- 
amond phase plate. Dynamical theory shows that bire- 
fringence occurs in perfect crystals when scattering at, or 
near, a Bragg reflection. Experimentally, the polarization 
ellipticity and handedness induced by the birefringence 
is determined by the thickness of the crystal and the de- 
viance from the Bragg angle. The efficiency of the 
phase plate is therefore limited by the beam divergence. 
Far from the Bragg condition the incident polarization 
is unrotated. As the Bragg condition is approached, a 
phase difference between pi and sigma channels relative 
to the crystal is induced and elliptically polarized light is 
produced. At a particular deviation a tt phase difference 
is induced and the phase plate behaves as a half wave 
plate. By maintaining this deviation from the Bragg con- 
dition and rotating on a x-circle about the beam, any 
orientation of linearly polarized light can be selected.— 

The polarization of the scattered beam was analyzed 
by a crystal selected and cut such that the beam Bragg 
diffracts with minimum absorption close to Brewster's 
angle of 45°; hence, only scattering x-rays with inci- 
dent polarization perpendicular to the analyzer scatter- 
ing plane. Therefore by rotating the analyzer crystal 
about the scattered beam by an angle a polarization 
selective measurement can be made. Poincare- Stokes pa- 
rameters were determined by fitting the integrated inten- 
sity of the analyzer rocking curve as a function of rj' to 
the following equation;^- 

/ = y [1 + PI cos 2r]' + P2 sin 2r]'] , (A.4) 

where rj' = 0° and 77' = 90° correspond to polarization 
perpendicular, a', and parallel, tt' to the scattering plane, 
respectively. 

By writing PI and P2 using Eqs. [2l and [A31 
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and 



IE, /I^i'2^(e'., e, z,, k', k)e^'»-"-J P + IE, %, k', k)e''>-J P 



P2 = — ■ — ■ — (A. 6) 



I 

where, as defined by Blume and Gibbs^ simulations of different magnetic moment directions and 

scattering factors can be made. By using a least-squares 
Zj = aj ui + hj U2 + Cj U3 (A. 7a) method, where a,, 6j, and Cj are free parameters, it is 

= U2 (A. 7b) possible to refine the magnetic moment direction from 



t^45o — [— sin 6> ui + U2 — cos 6 U3] (A.7d) 



_45o — — ^[sin6> ui + U2 + cos6> U3] (A.7e) 



= — sin^ ui — cos^ U3 (A. 7c) 

1 

e = sin(2x) sin^ ui + cos(2x) U2 — sin(2x) cos^ U3 

(A.7f) 

k' = cos6> ui — sin6> U3 (A.7g) 
k = cos6> ui + sin^ U3 (A.7h) 



a measurement of PI and P2 as a function of incident 
polarization. 
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